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Abstract 

We consider the gravity in the system consisting of the Bogomornyi-Prasad-Sommerfield (BPS) D3-brane embedded in the 
flat background geometry, produced by the solutions of the supergravity. The effective action for this system is represented by 
the sum of the Hilbert-Einstein and DBI actions. We derive the Wheeler-De Witt equation for this system and obtain analytical 
solutions in some special cases. We also calculate tunneling probability from Planckian size of D3-brane to the classical regime. 

PACS numbers: 11.25.Uv ; 98.80.Qc ; ll.25.-w 
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I. INTRODUCTION 

The discovery of D branes in the string theory have 
opened the way to the realization of the universe as a 
braneworld embedded in an ambient space. In the type 
IIA, IIB and I strings there are many D branes sustain- 
ing gauge and matter fields, at least one of them should 
include the Standard Model of particle physics. This D- 
brane will be correspond to the our universe. All the 
matter and gauge fields are due to the open string sec- 
tor. The strings from this sector have the end points 
constrained to move on the brane. The graviton and 
the dilaton are contained in the closed string sector and 
probe the extra dimensions. In the regions where the 
back-reaction of D-brane and stringy geometry can be 
ignored, the effective action for the worldvolume in the 
flat background is given by the DBI action coupled to the 
worldvolume gravity. Such models were used in cosmol- 
ogy in order to explain inflation from the string theory 

[!]• 

The DBI action describes the effects of virtual open 
strings at the tree level; it includes the effects of the 
background geometry and field strengths. But the non- 
linear form of the DBI action is inconvenient for quan- 
tization. However, introducing an intrinsic worldvolume 
metric one can write down the equivalent action [2] . From 
the other side one can express the DBI action as sum of 
constraints [3-5], where Hamiltonian for the D-brane is 
the constrain. In the first case, where DBI action is rep- 
resented by the intrinsic metric, one can consider this 
metric as a source of gravitation and the induced met- 
ric on the worldvolume as a field coupled to the grav- 
ity. This case will be considered. The physical sense 
of our study is that in the classical background of 10 
dimensional spacetime we consider 4 dimensional quan- 
tum subsystem. In the low-energetic approximation the 
classical backgrounds are given by the solutions of the 
supergravity. One of the best known and well investi- 
gated are AdSk. X Xio-k backgrounds (AdS/CFT cor- 
respondence). The other are backgrounds with warped 
metrices maintaining the four-dimensional Poincare sym- 
metry. The most interesting backgrounds are responsible 
for inflation and leading to the de Sitter vacuum in four 
dimensions [6]. 

In this paper we extend the applicability of the mod- 



els with the DBI action to the quantum regime. In this 
regime we have to take into account quantum effects as 
well as for the gravity and for the other fields from the 
string sector (or low-energetic approximation) . These ef- 
fects can be described by the Wheeler-De Witt equa- 
tion. Motivated by the above remarks we consider the 
Friedmann-Robertson- Walker (FRW) model coupled to 
the DBI action. In section 2 we derive a classical action 
for the considered system. This action is the base for 
section 3 where we obtain the Wheeler-De Witt equation 
and give solutions of this equation in particular cases. 
Section 4 is devoted to conclusions. 



II. DBI ACTION AND GRAVITY 

The well-known form of DBI action for a D3-brane is 
(modulo WZ terms): 
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(2.1) 

where ^ a p and B a p are pull-backs of a background metric 
qab and NS field Bab by an embedding field X, respec- 
tively. A strength F a p of U (1) gauge field is the intrinsic 
field on the worldvolume of D3-brane. As it is shown 
in [2] the above action can be expressed in the equiva- 
lent forms by an intrinsic auxiliary worldvolume tensor 
field h a p. We consider the simplest case when the fields 
F a p and B a p on the worldvolume vanish. Then the ac- 
tion (2.1) becomes the 4-dimensional Nabu-Goto, which 
takes the form: 
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where h a p is the intrinsic auxiliary metric on the world- 
volume and A is a constant. We promote the metric 
h a p to the dynamic degree of freedoms. Thus the met- 
ric h a p is a field which is responsible for gravitation on 
the worldvolume and the metric 7 on the worldvolume is 
considered as a field induced by the background. Hence 
the system consists of gravitation on the worldvolume 
coupled to the fields originated from the background. 
Problems of the backreaction in this consideration are 
negleced. So the action for this system is obtained from 
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(2.2) by adding the Einstein-Hilbert term: 



where 



d 4 x^- dct{h afj )R (h) - 

^ J d 4 xe-^-det(h a0 ) [h a ^ af3 - 2A] , (2.3) 

where nip = (8-kG)^ 1 . Next we apply the ADM con- 
struction for this system. In this construction the world- 
volume M is the product: R 1 x S 3 where S 3 is the 3- 
dimcnsional space-like slice of M . Then the metric h a p 
is determined by a shift vector N m and a lapse function 
N as follows: 



h 00 = -N 2 +h mn N m N n , h 



hmn-N , h rnn — h mnj 



where h mn is the intrinsic metric on S3. The matrix 
(/i Q/3 ) has the entries: 

h°° = -l/N 2 , h Qm = N m /N 2 , h mn = h mn -N m N n /N 2 . 
Using the above relations we get: 

h a(i la p = ~ (700 + N m N n lmn - 2N m l0m )+h mn lmn . 

The action (2.3) in the comoving coordinates (N m = 0) 
h m n — h mn and for the FRW metric : 



ds 2 = -N 2 dt 2 +a 2 (t) 
results in: 
S = 

2 
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(2.4) 
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where 

and 
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For our purpose we assume that S3 has O (4) symmetry 
and the metric h mn is the metric on S* 3 so k = +1. Hence 
the action reads: 



S = 6ir 2 m 2 P I dta 6 N 
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(2.6) 
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The function N(t) can be arbitrarily chosen by a redefi- 
nition of time thus we will use the gauge N = 1. 

This system we put in the flat background produced 
by N coincident BPS Dp-branes. The metric gMN, the 
dilaton <f> and RR-ficld C for this background have the 
form [7]: 

ds 2 w = g MN dX M dX N = B- x l 2 r] lxv dX> 1 dX v + Hl^dXjdX 1 



(fj,, v = 0, l...,p and I = p + 1, ...,9) . 



e 20 = H (3-p)/2^ 



C = (H- 1 - 1) dX° A ... A dX p , 



(2.7) 



(2.8) 



(2.9) 



where H p is the harmonic function of the transverse co- 



ordinates (AT/) to the world-volume: 

H p = l + ^, (r=(M J ) 1/2 ) (2-10) 



r f-p ■ 

and (rifi,,) = diag (-1, +1, +1). 

Here we consider backgrounds produced by £>p-branes 
with p > 3. We also use gauge freedom choice for an 
embedding field X in the form (so called static gauge): 

X(x) = (t,x 1 ,...,x 3 ,X 4 (t),...,X 9 (t)). (2.11) 

In this gauge the part of action related to RR field C 
vanishes. The induced metric in this embedding is: 
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where i = 4, ...,p. The integral on S 3 in the last term of 
(2.6) gives: 

J dne-+Tr{i) = 12tt 2 (l + 2 + 1 In 2^ H { p p - 5)/i . 
Thus the considered action (2.6) takes the form: 
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where L is given by the formula: 
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and 7 = 1 + 2 + i ln2 = 3.3466. In the directions spanned 
by Xi and Xj we introduce spherical coordinates: 

Xi = gfi, Xi = rhi, 

where fif 1 — hjh 1 = 1. In the generic case the consid- 
ered system (the probed £>3-brane) has a non-trivial an- 
gular momentum in the transverse directions (Xi) and 
non-trivial angular momentum in the directions trans- 
verse (Xi) to T>3-branc but longitudinal to the back- 
ground branes. For the simplicity, we consider the non- 
rotating case: fi = hi = 0. Thus, the cq. (2.13) takes 
the following form: 



In the case when r is fixed, the equations for a are 
obtained from the Hamiltonian constraint and from the 
equation of motion for (2.17). These equations are given 
by: 
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Then the three dynamic fields (a, g, r) span minisuper- 
space, and the action is: 
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If one compares these equations to the standard Fried- 
mann equations for a perfect fluid with an energy density 
p, a pressure p and a curvature parameter k: 
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where a potential U is given by: 

U (a, r) = a 3 AT 3 (h^/ 4 - 2A)-a (m 2 P + AT^H^^) 

(2.15) 

The metric Gs$ on the minisuperspace is taken from 
(2.14): 
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(2.23) 
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AT 3 a 3 i4 p - 1)/4 



(2-16) 

The equation of motion for g gives the following relation: 



/ then one can notice that the energy density p is given by: 
p = 3AT 3 (h^I 4 - 2A) + A J 2 ^ 5 "f)/ 4 , (2.24) 
the pressure p is expressed as follows: 

p = -3AT 3 (H^l 4 - 2A) + A j2^5- P )/4 (2 25) 
and the curvature parameter fc is given by: 
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In this way one obtains: 

where J is a constant. Thus eq. (2.14) is: 
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(2.17) 



The first term on r.h.s in (2.24) can be interpreted as a 
cosmological constant A : 

A = 3AT 3 (h^-V' 4 - 2A) . (2.27) 

The state equation for this perfect fluid has the form: 



and a potential U is given by: 

U (a, r) = a 3 AT 3 (h^' 4 - 2A 
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a ( m 2 -|- AT }j(p- 5 )/ 4 \ _|_ J 2 a _3 ij( 5_p )/ 4 (2 18) H cnce for a ^ oo one obtains the following state equa- 
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Now, the mini-superspace is reduced and spanned by a 
and r. The metric Gs$ on this space is: 



P 

w = — 



-1. 



(Gs$) — 
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Thus, for big a the worldvolumc of D3-brane will be dom- 
inated by the perfect fluid with the negative pressure and 
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the observer fixed to this worldvolume will sec an accel- 
erated expansion. 

Let us notice that for p = 5 the action (2.17) is: 



= 6tt 2 / 

JR1 



dt 



-m 2 P aa + AT 3 a 3 H 5 r 2 + 
a (m 2 P + AT 3 a) - J 2 a" 3 - a 3 AT 3 (1 - 2A)] , (2.30) 
where 
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If one replaces the field r by the following field ip: 



Lp = x/Ng s + r 2 - Ng s In 



Ng a + V N 9s + r 2 



(2.31) 



then the action takes the form: 



6tt z 



dt 
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-mpaa + AT 3 a 3 p 2 + 



a (m 2 P + AT 37 ) - J 2 a~ 3 - a 3 AT 3 (1 - 2A)] . (2.32) 

Thus for the background produced by 5-branes the ac- 
tion for 3-brane is reduced to the free scalar field ip 
coupled to the scale factor a with the potential v — 
-a (m 2 P + AT37) + J 2 a- 3 + a 3 AT 3 (1 - 2A). For J = 
the action (2.32) is the action for the homogeneous and 
isotropic metric g coupled to the free scalar field with the 
cosmological constant AT 3 (1 — 2A). One can see that 
the cosmological constant A is the function of intrinsic 
parameters describing 3-brane, namely the strength T 3 
and the constant A. Thus by a suitable choice of A one 
can obtain the real value (observed) of the cosmological 
constant. For p ^ 5 the cosmological constant is given 
by (2.27). 



III. THE WHEELER-DE WITT EQUATION 

In the mini-superspace spanned by the fields Q = 
(Qi, Qn) the general form of the Wheeler-De Witt 
(WD) equation has the form: 



h 2 1 



:d @ (7^GG en d n ) +U(Q) 



* (Q) = 0, 



2 y/^G~ 

(3.1) 

where Gen is a metric on the mini-superspace. In our 
case the metric is given by (2.19) and Q = (a, r). In the 
equation (3.1) one has to take into account the factor 
ordering between the conjugate variables. The kinetic 
terms of the Hamiltonian for a and r are obtained from 
(2.17) and can be expressed in the equivalent forms on 
the classical level: 

^ = a-^ +k+ ^P a aiP a a k , 



TT- [{p-Vj / i+U+w] p ETttp VfW 

n p r rn p r r n, p , 



AT 3fl 3^ p - 1)/4 AT 3 a : 

where j, k, u, w are arbitrary numbers. On the quantum 
level they are non-equivalent and depend on the commu- 
tation relations between the conjugate variables. This 
ambiguity has to be respected in the quantum Hamilto- 
nian. In the canonical quantization the ordering ambigu- 
ity in the expression / (x) P 2 for the conjugated variables 
(x,P x ) is resolved as follows: 



f(x)P 2 ^-h 2 (l + A + B)fx 

2B f d B f" 



d 2 _a 



A + B f dx 1 + A + Bf 



where A, B are constants. Hence the differential part of 
the Hamiltonian operator in our problem reads: 
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4H n 
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where 7, 6, /i, v represent the factor ordering ambiguity 
and prime denotes derivative with respect to r. Thus the 
WD equation is: 
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(3.2) 



where an effective potential U e f / (a, r) has the form: 

U eff (a, r) = a 4 AT 3 (h^'± - 2 a) - a 2 (m 2 P + AT 37 ^ p - 5)/4 ) + 
J 2 a- 2 Hf-rtl\ (3.2a) 

At the moment we assume that r is fixed and is consid- 
ered as a parameter. Thus the potential U e ff and the 
wave function ^ have fixed values for r so the eq. (3.2) 
results in: 



da 2 
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2m%k o 



2m 2 p \ 4 



* = 0, 



(3.3) 

where k, A are given by (2.26-2.27) and r\ = —C — 5 + vD 
with C and D given by: 

_ 2m 2 p 2 ( 5 _p)/4 

- J n v 
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D = 



m 2 P (p - 1) 
4AT 3 ^ p+7)/4 



{(p + Z)H' p -AH p H'^. 



Then for the fixed value of r the equation (3.7) becomes 
the Bessel equation: 



The last coefficient D has the form: 



D(r) = - 



mp (p — 1) 
4AT 3 i4 P+7)/4 



9(r) 



(3.4) 



d 2 F 1 dF 
dz 2 z dz 
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^0 = 0. 



(3.10) 



The solution of (3.10) is given by: 



and: 



_ 4(7-p)(8-p)Ng 3+p ,1^5- 



r 8-p 



L 4(8-p) r r 8 -P J 



(3-5) 

The eq. (3.3) can be transformed to the following form: 



F (z) = EJ n (z) + FY n (z) , 

where J n and Y n are the Bessel functions of the first kind 
of the order n. Hence the wave function \& obtained from 
(3.10) has the form: 



d 2 F 4s + 1 - 7 dF 



dz 2 



+ 



2z dz 



+ (l- — +mz)F = 0, (3.6) *(a;r ) -a (1+7)/2 



where a variable z is related to a as follows: 



z = —la 



2 m 2 P I k 



H V 2 



and the wave function \I> is expressed by a function F in 
the form: 



*(z) = z s F(z). 
The coefficients m and n 2 are given by: 



(3.7) 



m = —i ■ 



HX 

3mpk V fc' 



^ 2 = t 7 -2a(2s-7-l) 



If we choose s = (1 + 7) /4 in (3.6), then this equation 
becomes: 



d 2 F 

dz 2 



where: 



IdF / n~ 

--r + 1 1- — 



n 2 =r ? /4+( 7 +l) 2 /16 







(3i 



(3.9) 



This equation is similar to the Schrodingcr equation con- 
sidered in the Stark effect [10]. Because this equation 
has non analytical solutions the wave function for the 
Stark effect is given only in a perturbation series in m. 
First we consider the case when m = with a solution 
F (z). The solution F (z) corresponds to vanishing of 
the cosmological constant A with the mirage perfect fluid 
on the worldvolume. This fluid has the energy density 
p = 3 J 2 Hp 5 p ^ 4 /a 6 and the state equation w = +1 (sec 
eqs.(2.24) and (2.28)). Under our assumption that r is 
the parameter we are able to fix it in such a way that m 
is equal to zero. Thus one can relate A to the position 
r : 



2A = #(p- 5 )/ 4 (r ) 
and the parameter k takes the form: 



(3.9) 



k= 1 + 



EIr, 



l pi 



FK n ( 



(3.12) 

where E and F are constants, /„ (z) and K n (z) are the 
modified Bessel functions of the first kind and the sec- 
ond kind. We introduced the Planck length lpi = H/mpi. 
This wave function must satisfy given boundary condi- 
tions. It has to be regular everywhere, so using expansion 
of /„ and K n near zero we obtain: 



1 + 



7 > ^477+ (7+ l) 2 



It means that 4 77 < 0, so in explicit form we get: 



2™^ 



J 2 - 5 + vD (r ) < 0, (3.13) 



where D (ro) is given by the equation (3.4). From this 
relation one can, in principle, determine the allowed po- 
sitions of D3-brane in the background for which the wave 
function of the D3-brane is given by (3.12). For z — > ioo 
the asymptotics for I n and K n are following: 



K n (z) 



'2t:z 



exp (z) , 



exp (—z) . 



Thus the wave function (3.12) for big a reads: 



* (a; r ) ^ lpi 



\n 2 k 



1/4 



,(7-l)/2 



I 2 
L Pl 



1 PI 



eexp ( + — v / H 2 k/2 
(3.14) 



The effects of the non-trivial D-brane background are in- 
cluded in k and in the order n of the Bessel functions. 
The regularity of the wave function near zero puts the 
constraint (3.13) on the positions of the D3-brane in the 
background. The predictions following from this wave 
function are well-known and broadly discussed with re- 
spect to the decay of a false vacuum [9] . 
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The second exact solution of (3.6) is obtained for s = 
(7 — 1) /4 and n 2 = 0. In this case (3.6) takes the form: 



d 2 F 

dz 2 



+ (1 + mz) F = 0. 



(3.15) 



It is then transformed to the Airy equation: 
dt 2 m 2 



where t = 1 + mz and F (z) = 4> (1 + mz). The solutions 
$ of the Airy equation depend on the sign of m 2 . For 
m 2 > and for m 2 < the solutions are given by: 



*+ (*) 



1/3 



3 m 



i 3 / 2 + B' J 



1/3 



3 TO 



7 i 3 /< 



(3.17) 



•-(')- ^ 



^-1/3 



3/2 



B7 



ol I " / 1 — -1/3 1 ol 

3|m| / ' \3|m 

' (3.18) 

respectively. Since to 2 = — ( J (7^) < the so- 
lution will be given by the second relation. Finally the 
wave function ^ has the form: 



.3/: 



* (a; r ) 




Bh/3 



l A Pl \ 



(3.19) 



and s = (7 — 1) /4. This wave function is valid only for 

2 -0: 



n 



»7+(7-l)(7 + 3)/4 = 0. 



(3.20) 



The above condition fixed the position ro of D3-brane in 
the background because r\ is given by H p (see eqs. (3.4) 
and (3.5)). The wave function (3.19) has the oscillating 
character for 3h 2 k/ (AZpj) < a 2 and for all a > is regu- 
lar. The (3.19) has the form of the Hartlc-Hawking wave 
function for empty universe with the cosmological con- 
stant [11]. It is as should be since the condition (3.20) 
implies that the worldvolume is empty with the cosmo- 
logical constant A. 

Thus we have obtained two particular solutions of the 
cq. (3.6) which correspond to the fixed positions of the 
D3-brane in the background. The first solution (3.12) 
is obtained for vanishing of the cosmological constant A 
on the worldvolume which is filled by the mirage fluid 
with the energy density 3 J 2 Hp 5 p ^ 4 /a 6 . This solution is 
valid only in the region where the condition (3.13) holds. 
The second solution (3.19) is obtained for the position of 
D3-brane fixed by (3.20) and corresponds to the empty 
worldvolume with the cosmological constant A. 



Next we consider the tunneling probability in the po- 
tential (3.2a) for a fixed position ro of D3-brane with 
to, n 7^ 0. In order to do this we choose s = (7 — 1) /4 in 
(3.6). Thus we get the one dimensional equation: 



d 2 F 

dz 2 



2 



0, 



(3.16) where: 



n 2 = 77/4 +(7-1) (7 + 3) /16. 



(3.21) 



(3.22) 



This equation can be considered as a Schrodinger equa- 
tion in a potential V (z) with a zero eigenvalue: 



d 2 F 



where V (z) is given by: 

V(z) = -1 



mz. 



(3.23) 



(3.24) 



For the real variable y = Jh 2 k/2 < so that z = iy 

l pi 

the potential V takes the form 



V(y) 



-1 



(. 



f +xn \h 2 k 



3/2 



(3.25) 



where x — Zp ; A/ (6n 2 ) and In is the Planck length. This 
potential is sketched below: 



The variable y has physical meaning only for y < 0. The 
potential V has a maximum for 



1/3 



ffk 



with the value: 



-1 - 3n- 



(I) 



2/3 



1/2 



2 



(3.26) 
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The scale factor a 2 , which corresponds to j/ max is: 



I 2 
1 pi 



1/3 



and is grater than Ipi for x < 2. In the classical regime, if 
an energy E of the system is less then F max , we never pass 
from the region where a < a to the region where a > a . 
It means that the D3-brane collapse to the zero size. But 
in the quantum regime exists non zero tunneling proba- 
bility to the region where a > a which corresponds to 
expansion of D3-brane worldvolume. This probability T 
in the quasi-classical approximation is given by: 



T ~ exp 



2 f y3 n 



E + V (y)dy 



(3.27) 



where the integration limits 2/3, 2/2 are given by the solu- 
tion of the equation: 

E - V (y) = 0. 

The roots 2/3, 2/2 are negative and ordered as follows: 2/3 < 
2/2 < 0. Hence for y E (2/3,2/2) and E < V max < we 
obtain that: 

E - V (y) < 0. 



3/4 



Thus, an integral / = J^ 3 —E + V (y)dy is real and has 
the form: 

I=nV ~ x (ik) JV2 

a function / is given by a cubic polynomial: 

/ (v) = v 3 + °y 2 - w, 



Jm V 



(2/), 



(3.28) 



(3.29) 

where a = -w (E + 1) / (n 2 ) and w = (fi 2 fc/2) 3/2 fx. 
Now we have to determine the integration limits 2/3and 2/2 
which depend on E. We make the reasonable assumption 
(in the quasi-classical approximation) that the minimal 
value of the scale factor a on the worldvolume is given by 
the Planck length Ipi = h/mpi. This minimal value cor- 
responds to y pi = —hy/k/2. At this point the potential 
V has the value: 



V(ypi) 



-1 - 2n 



2 1 + x 

h 2 k ' 



and Vmax > V (ypi). The coefficient a for E 
takes the form: 



(3.30) 
V(ypi) 



a = 



h 2 k 



1 + 



1 



(3.31) 



We evaluate Y near the energy E = V (ypi). Thus: 



y/-E + V(y) = y/-V (yp,) + V(y) = 

( 2 \ 3/4 1 
= n ^\Wk) y^- 



(3.32) 



It means that the cubic polynomial / (y) has a root for 
y = ypi < and the bottom integration limit is 2/2 = ypi- 
The other roots are given by: 



yi 



2x 



h 2 k 



(-i + >/r+4s) 



(3.33) 



1 h 2 k 



V3 = '2i\~ (l + VYTI ^ (3 ' 34) 

and they are ordered in the following way 2/3 < ypi < 
and 2/1 > 0. The upper integration limit is 2/3 and the 
polynomial / = (y - y 3 ) (y - y Pl ) (y - 2/1) is positive for 
V <= {V3,ypi)- The barrier width b which is equal to 

ypi - 2/3: 



b = d —{- 1 + — 2~ x — 



(3.35) 



is greater then zero for x E (0, 2) and vanishes for x = 2. 

The integral (3.28) is the elliptic integral and is ex- 
pressed as follows: 



\ 3/4 



where: 



-l 



dy 



(3.36) 



(3.37) 



and s = —1,1. Making the standard substitution y = 
2/3 + (ypi — 2/3) sm2 m the integrals I 8 we obtain: 

h = - 7 ^= K( K ) + 2Vy7=wE(K), (3.38) 
V2/i - 2/3 



2/3\/2/l - 2/3 



n ( c , «) , 



(3.39) 



where K (k) , _E (k) and n (c, k) are the first kind, the 
second kind and the third kind complete elliptic integrals, 
respectively (see Appendix). The modulus n and c are 
given by: 



ypi - 2/3 1 f 1 . 1 - 2.x 
K = = ^ H ; 1 

2/1-2/3 2 V VI + 4xy 



2/P( - 2/3 1 



(3.40) 



2/3 



- (3 - VT+ Ax) > 0. (3.41) 



They have the range: < n < 1 and < c < 1 for 
x E (0, 2). Thus the integral I takes the form: 



ny/x 



( — 

\H 2 k 



3/4 



2o-2/i 20-^/^3! 

=A (k) H 7- E (k) ■ 



3\/2/3 



2w 



2/3V2/1 - 2/3 



n(c, k) 



(3.42) 
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where 1/31 — y\ — 1/3 > and j/32 = jj2^V3 > 0. Collecting 
the above formulas we get the function I in the variables 
x, and n: 

(1 + 4x) i/4 

n 1 + X [AVTT^E (k) + (1 - >/l + 4£) X («)] . 

6a; (1 + 4a:) 7 

(3.43) 

This function tends to 00 if x — > and is decreasing 
function of x with the minimum for a: = 2 

I(2,n) = ^=. (3.44) 

For the value of a; = 2 the barrier width 6 vanishes and 
n 2 = Zp;A/12. The maximum of the tunneling probabil- 
ity T is given by: 

T max = exp(--/ (2, n)) = cxp ^"^J • (3-45) 

Our result agrees with [12] in the sense that maximum of 
the probability is achieved in the point where the barrier 
potential vanishes. Moreover n is the function of r and 
parameters related to the ordering (see eqs. (3.22), (3.4) 
and (3.5)). Thus r max takes the maximum for n — with 
the cosmological constant A = on the world-volume. 
This case corresponds to the position r max for which: 

n (r max ) = 0. 

Thus r max gives the boundary value of the condition 
(3.13). In this position the D3-brane has the wave func- 
tion given by (3.12). 

IV. CONCLUSIONS 

We have considered the gravity on the worldvolume 
of D3-brane embedded in the flat background produced 
by N p-branes. This system is described by the sum 
of the Hilbert-Einstein and DBI actions. Although the 
DBI action is non-linear we have used equivalent form for 
it with an auxiliary metric. This metric has been pro- 
moted to the dynamic field on the worldvolume. Thus 
the worldvolume gravity is described by this metric. In 
the classical regime we have obtained potential (2.18). 
Comparison of the classical equations of motion for this 
potential (for fixed r) to the Friedmann equations results 
that (2.20) and (2.21) can be interpreted as the energy 
density and the pressure of the mirage perfect fluid. For 
the big value of the scale factor a we get the state equa- 
tion of this perfect fluid in the form (2.27). Thus one 
can conclude that the expansion of the worldvolume is 
accelerated. 

In the next part of this paper we have derived the 
Wheeler-de Witt equation for this system and solved it 



for the special cases corresponding to the fixed positions 
of D3-brane in the background. In the first case the 
wave function has the form of the instanton tunneling 
and in the second we get the Hartle-Hawking function. 
The obtained wave functions depend on the ordering of 
the conjugated variables. In general, the wave function 
depends on a and r but we can not find it in the explicit 
form. Even for fixed r we get an equation which has not 
explicit solutions and is similar to the equation consid- 
ered in the Stark effect. The perturbation method used 
in the Stark effect is inappropriate because we do not 
know which parameter is small in eq. (3.8). We made 
the assumption that the minimal value of the scale factor 
a of the worldvolume is given by the Planck length and 
for this value we evaluated tunneling probability in the 
potential V. This probability is finite and depends on A 
and r. The maximum of the probability is for x = 2 and 
the cosmological constant is given by A = I2n 2 lpf. 

The considered model has the common features with a 
hydrogen atom in an external electric field. In our model 
the external field is represented by the background pro- 
duced by Dp-branes and the atom is represented by D3- 
brane. The tunneling from Planckian size of D3-brane 
to classical regime corresponds to the ionization of the 
hydrogen atom. 

The above model is very simple and not quite realistic 
for several reasons. The first one is that we used FRW 
ansatz for the metric. Validity of this metric was ex- 
tended to the pre-inflation epoch, when quantum gravity 
has dominating effects. In this regime the universe was 
not necessary isotropic and homogenous. The better ap- 
proximation for the metric is IX Bianchi type. This case 
were considered in [13] for the empty universe with a cos- 
mological constant. The second reason is that the stringy 
effects enter by the DBI action. This action is valid only 
on the tree level of the string theory and does not take 
into account backreactions. The other reasons are re- 
lated to the special backgrounds which are classical so- 
lutions of the supergravity equations. The more general 
approach is if one consider the Wheeler-De Witt equa- 
tion for the low-energetic approximation of string theory 
by supergravity and next compactificate this equation to 
4-dimensions taking into account stringy corrections. 

In the classical regime compactification to 4- 
dimensions for the Bianchi type-I cosmology in the pres- 
ence the gravity, the dilaton, and the antisymmetric ten- 
sor field of the second rank, coupling to the gauge field 
strength living on the D3-brane was considered in [14]. It 
would be interesting to consider the quantum regime of 
this system, described by the Wheeler-DeWitt equation, 
and compare the results with the results of the model con- 
sidered in this paper. We shall investigate this problem 
in the future. 
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V. APPENDIX 

The complete elliptic integrals first kind K, second 
kind E and third kind II are defined as follows [15]: 

tt/2 



Jo y 1 — K SI 



,tt/2 

E(k)= / - /,-siir 

JO 



II(c,k) 



tv/2 



dtp 



o (l — csin 2 </>) y/l — 



K sin 



For c — k the elliptic integral II is expressed by E: 

E{k) 



n (k, k) 



1 - K 
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